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Abstract. We review recent results on the charged pion - charged hadron correlation in p + p 
and d + Au collisions as measured by the PHENIX Collaboration. Properties of di-jet system, 
such as the jet shape, associated hadron yield per trigger pion, and the underlying event are 
extracted statistically from the tt^ — correlation function in and A77. For jet triggered 
with high pt pions (pr > 5 GeV/c), no apparent differences in the jet properties are seen 
^~») ' ■ between p + p and d + Au. 
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1. Introduction 

^ The technique of two particle correlation in relative azimuth (A(f>) and pseudorapidity (A77) is 

an useful tool to access the (di-)jet properties in heavy-ion collisions. Comparing with the 
traditional full jet reconstruction method, the two particle correlation method is relatively 
insensitive to the level of the underly event, thus can probe soft jets (<5 GeV/c); combining 
with event mixing technique, it can also be used for detectors with limited acceptance. 

To leading order in QCD, high px jets are produced back-to-back in azimuth. This back-to- 
back correlation, however, is smeared by the fragmentation process and the initial and final state 
radiation, which lead to a typical di-hadron correlation function in A<p as shown schematically 
in Figure^ The associated hadron yield per trigger (conditional yield or CY) can be 
parameterized by a constant plus a double gauss function, 
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In this analysis, everything about the (di-)jet is extracted from this parameterization. The 
peaks in the same side (Acj) = 0) and the away side (A(p = tt) represent the intra-jet and di- 
jet correlation, respectively. The widths of the peaks are controlled by the jet fragmentation 

momentum Jt and the parton transverse momentum kx o sarne cx Jt v , <J ^t way ~ a 1ame ^ 

kxyi where the subscript "y" represent the ID projection in transverse plane; The integrals of 
the peaks, Ns and Na give the total number of hadrons associated with the trigger hadrons in 
the same side and the away side; The pedestal level beneath the jet structure, B, represents 
contributions from the underlying event. 

Thepresence of the medium can modify the di-hadron correlation. Recent results from 
RHIC H Q, and SPS Q indicate a complicated modification of the jet structure in both 
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Figure 1. Cartoon of the two particle Ac/) correlation. The yield of hadron per trigger 
(Conditional Yield) has a di-jet part and a part corresponding to the underlying event. 



the same side and the away side. The same side jet become elongated in A77 but is relatively 
unmodified in A(f> direction. Meanwhile, the away side correlation become much broader and 
suppressed, which indicates strong interaction of the jets with the dense medium. In order to 
achieve quantitative understandings of the modifications of jets in the medium, one has to obtain 
an accurate baseline measurements of jet correlation in p + p and p + A collisions. Di-hadron 
correlation in p + p collisions probes basic QCD effects such as jet fragmentation, initial and 
final state radiations; while correlation in p + A collisions gives handle on the various initial state 
effects such as shadowing and jet broadening in cold nuclear medium. 

We focus the physics discussions on three aspects of the tt^ — h correlations in Figure^ 
jet shape, jet yield and the underlying event. But before doing that, we briefly discuss the 
identification of high pr charged pion and the techniques used to extract the jet properties. 



2. Analysis 

2.1. tt^ identification 

PHENIX identifies high momentum charged pions with the RICH and EMCal detectors. 
Charged particles with velocities above the Cherenkov threshold of 7^ = 35 (CO2 radiator) 
emit Cherenkov photons, which are detected by photo-multiplier tubes (PMTs) in the RICH 0. 
This threshold corresponds to 18 MeV/c for electrons, 3.5 GeV/c for muons and 4.9 GeV/c for 
charged pions. In a previous PHENIX publication ||, we have shown that charged particles with 
reconstructed pt above 4.9 GeV/c, which have an associated hit in the RICH, are dominantly 
charged pions and background electrons from photon conversions. The efficiency for detecting 
charged pions rises quickly past 4.9 GeV/c, reaching an efficiency of > 90% at pt > 6 GeV/c. 

To reject the conversion backgrounds in the pion candidates, the shower information at the 
EMC is used. Since most of the background electrons are genuine low pr, they can be rejected 
by requiring a large shower energy in the EMC. In this analysis, a momentum dependent energy 
cut at EMC is applied: E > 0.3 + O.lbpx- Additional electron rejection comes from the x 2 
variable, 
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where ^? ieas is the energy measured at tower i, Ef red is the predicted energy for an 
electromagnetic particle of total energy Yli Ef 1 ^ and o~i is the standard deviation for Ef Ted . 



Both Ef ved and Oi are obtained from the electron test beam data. EM shower is more compact 
than hadronic shower, thus has a smaller \ 2 value. The x 2 value is then mapped to the 
probability (prob) for a shower being an EM shower, prob ranges from to 1, with a flat 
distribution expected for EM showers and a distribution peaked around for hadronic showers. 
Figure. [2^ shows the normalized prob distribution for the pion candidates and electrons. A cut 
of prob < 0.2 selects pions with an efficiency of > 80%. Since we are interested in per-triggered 
yield, the detailed knowledge of the pion efficiency is not necessary. The raw pion spectra for 
requiring only RICH cut and both RICH and EMCal cuts are shown in Figure. HJa, the difference 
between the two is mostly due to electron background. The sample of charged pion used in the 
correlation analysis is from 5 to 16 GeV/c, with an purity better than 95%. 
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Figure 2. a) The probability distribution for charged pion candidates (solid circles) and 
electrons (open circles), where the integrals of the histograms have been normalized to one. 
b) Raw pion spectra with matching hit at RICH only (solid circles) and with both RICH and 
EMCal cuts (open circles). 



2.2. Extracting jet properties 

The correlation functions are generally defined as 

iVcor(A<A) ■ . Ar <(\\ ^or (At/) 

C A< ^ = Tt 7X^\ m azimuth and C ( Ar ?) = 7a~T ( 3 

iVmix (A0) iV mix (At?) 

in pseudorapidity. A^ COI and iV m i x represent the same-event pair distribution and mixed-event 
pair distributions, respectively. The mixing is done by pairing trigger tt^ with charged hadrons 
from events having similar collision vertex and centrality as the -k^ 1 . Ref.[3] has shown that the 
correlation function and the conditional yield are related to each other by just a constant, 
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where N® T - and JV t rig are the true and detected number of triggers respectively, and e is the 
average single particle efficiency for the associated particles in 2tt in azimuth and ±0.35 in 
pseudo-rapidity. accounts for the loss of jet pairs outside PHENIX pair acceptance of 

\At]\ < 0.7. 

The two gauss functions in Eq^ describe the A(j) distribution of the jet signal. The jet signal 
can also be presented in any other pair variables, such as Arj, trigger px (pr,trig); associated 
hadron p T ( pr.assoc), Pout = PT,assoc sin A(p, x E = PT - a °t° cCOsA<? \ di-hadron mass and di-hadron 

¥ 1 .trig 



Pt- For every pair variable, we use a statistical weighting method to account for the acceptance 
correction. According to Eq0J each pair on average needs a Acj) dependent correction factor, 
w(A<f>), 

S / dA</,JV mi x(A</>) 

When this factor is used as the weight in filling the xe histograms for both real and mixed pairs, 
we obtain 
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for the same-event pair distribution. Thus according to EqQ X E distribution for jet pairs equals 
to 

1 r7/V jet 

jyo. dx = X) *(«*>(Afl S(x E )w(Ad>) . (7) 
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where 

c = BR Av 2tt 

N trig e J dA(f>N mix (A<f>) 

Replacing xe with any pair variables, we obtain other jet pair distributions. However, the 
integral of the jet yield should be conserved independent of the pair variable used, i.e.: 

dA /At = I ix /J^ = j dPTm J^L = I (9) 

dA(p J dx E J CtPT.assoc J dp out 



3. Results 

3.1. Jet shape 

In the following discussion, the trigger -k^ px is always from 5 to 10 GeV/c, unless specified 
otherwise. FigureOl shows the tt^ — Acf) distribution from p + p and d + Au collisions for 
several range of pt, assoc- The widths decrease with increasing pt, assoc, which is consistent with 
narrowing of the jet cone for larger pt, assoc- It is interesting to notice that a large fraction of all 
hadrons in the event are associated with the trigger, thus are originated from the hard-scattered 
partons. Even for pT,assoc as low as 0.4 — 1 GeV/c, about 51% hadron yield in p + p (27% in d 
+ Au) comes from the jet fragmentation. 

Using the event mixing technique, we also measure the jet shape in rj. For the single 
acceptance of \r]\ < 0.35, the pair acceptance in pseudorapidity is limited to be \Ar]\ < 0.7. 
Figured shows the same 6V6nt cincl mixed, event At] distribution for 1 <C PT,assoc 

< 2 GeV/c, 

where a cut of \A(f)\ < 1 is used to select only same side jet pairs. The mixed-event pair 
distribution is not a perfect triangle due to a gap around rj = in PHENIX central arm 
detectors. The ratio of the two distributions gives the jet shape in Arj. It is shown in FigureEJ) 
and compared with the jet shape in A(p. There is no significant difference between the two 
and the extracted widths are consistent in both directions. We extend this comparison to other 
associated hadron pt ranges and summarize the results in FigureEl The overall agreement 
between the jet widths in A77 and Acj) is pretty good, except at small PT.assoo where the width 
in Arj is systematically lower than that in A(f>. The fact that this discrepancy exist in both p+p 



< 
-a 



1.5 
1 

0.5 



0.2 

-e- 
< 

"D 

= 0.1 



+d + Au 

*P + P 



PHENIX Preliminary 



0.4<p_. <1.0 GeV/c 

n T,assoc 



CCD- 



2-0<p T „„„„<3.0 GeV/c 



i 

r • ? -^0° A. i 



2 
A4>(rad) 




2 4 

A^(rad) 



Figure 3. Corrected conditional pair distributions for p+p and minimum bias d+Au collisions. 
The trigger -k^ are correlated with hadrons with pr jassoc 0.4—1.0 GeV/c, 1.0—2.0 GeV/c, 2.0—3.0 
GeV/c and 3.0 — 5.0 GeV/c (from top to bottom and left to right). 




Figure 4. a) The same-event and mixed-event pair distribution in Ar/, b) the correlation 
function in Arj (open boxes) and A(j> (filled circles). 



and d + Au collisions indicates that this deviation is likely due to the systematics of the fitting 
in a limited Ar/ range rather than any real physics effect in d + Au. 

From the measured jet width, one can calculate the rms value of zkx y [2], ( z kTy)nMS = 

(^z 2 k Ty ^, for both p+p and d + Au. The resulting (zkTy) rm s is plotted as function of trigger 

Pt in FigureEl It looks quite similar between p + p and d + Au. The difference of (zk r py)RMS^ 
averaged over px, is {z 2 k Ty ) dA u ~ {z 2 k Ty ) pp = 0.64 ± 0.78 ± 0.42 (GeV/c) 2 , which is consistent 
with 0. According to various theoretical estimations jjj, the typical contribution to (k Ty ) from 
multiple scattering is 1 (GeV/c) 2 in central d + Au collisions, while the contribution from initial 
and final radiation is much larger (around 8 (GeV/c) 2 ) [lfj. The small multiple scattering 
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Figure 5. The comparison of jet width as function of pr,assoc m Acj) (solid circles) and Ar] (open 
boxes) from 7r ± — correlation, a) results for d + Au. b) results for p + p. 



contribution might have been washed out by the much larger contributions from initial and final 
radiation, which explains the lack of difference between the two systems. 
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Figure 6. The comparison of (^trig^T y )RMS 
values between d + Au (filled circles) and p + p 
(open circles) as function of PT,trig- 
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An alternative but more direct way in studying the away side kx broadening is through the 
p ou t distribution. For small angles, p ou t has simple relation to jx and &;t[3] : 



(Pout ,samei 



(jT,y) + x E(jT,y) 



(Pout,away) ~ 0t,j/) + x E(jT,y) + 2x%(z 2 k Ty ) 

Thus the difference of same side and away side p ou t distributions directly reflects the contribution 
from kx- 
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Figure. [TJi shows the same side and away side p ou t distribution. There is a significant difference 
between the two, which reflects the contribution from kx- Both distributions have a gauss shape 
at small p ou t followed by an excess at large p ou t- The gauss part presumably is due to the jet 
fragmentation (in both the same and away side) and intrinsic kx (away side only), while the 
excess is evidence for hard radiation contribution of the outgoing partons. Since the away side 



Pout carries information about kj<, we compare between p + p and d + Au to see whether there 
is hint of additional kx broadening in d + Au. The comparisons of away side p ou t distributions 
are shown in FigureEb> no apparent differences are observed, consistent with the observations 
that (zkT y )RMS are similar between p + p and d + Au. 
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Figure 7. a) Same side and away side p on t distributions in d + Au collisions, b) The away 
side p ou t distributions compared between p + p and d + Au. 



3.2. Jet yields 

The same side and away side px distributions of the charged hadrons associated with jets are 
plotted in Fi^Hl comparing between p+p and d + Au collisions. The same side yield is related to 
the di-hadron fragmentation, since both particles comes from the same jet, while the away side 
yield depends on two independent fragmentation functions: one parton fragments to produce 
the trigger, while the second parton produces the associated hadron. No apparent differences are 
seen between p+p and d + Au; this observation is in contradiction to some recombination model 
prediction [TJ|, in which a significant difference is expected due to shower-thermal contribution. 
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Figure 8. Jet pair distribution as function of PT,assoc for same side (right panel) and away side 
(left panel) in p + p and d + Au. 



Di-hadron correlation also gives the xe distribution 
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di-jet pt imbalance is ignored and z m 
xe distribution is closely related to the fragmentation function D(z) 



, where x E = ZassocAtrig- When 
•C ^trigj ztrig varies very slowly with z a ssoc- Hence the 
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Figure El shows the measured xe distribution between p + p and d + Au. Again, no difference 
is seen between the two collision systems in both the same side and away side. 
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Figure 9. xe distributions as function of PT,assoc for same side (right panel) and away side 
(left panel) in p + p and d + Au. 



In e + e~ or p + p collisions, the fragmentation functions D(z) are known to approximately 
scale, i.e. are independent of jet energy. To check whether this is still true in d + Au collisions, 
we plot in Figure. EH the conditional yields as function of trigger px in different ranges of xe 
for both p + p and d + Au. The amount of variation is within + 25% for pt from 5 — 12 GeV/c 
with very little difference between the two systems. So we conclude that the evolution of the jet 
fragmentation as function of jet energy is very similar between p + p and d + Au. 

3.3. Comments on underlying event 

Events triggered by high px hadrons not only contain particles originated from the two hard- 
scattered partons, but also those come from soft multiple interaction and the beam remnants. 
Underlying event in p + p and d + Au collisions refers to all hadrons except those from the 
two outgoing hard-scattered partons, which includes contributions from the beam remnants and 
initial and final state radiation ^3] • The physics of the underlying event is poorly known due to 
its non-perturbative nature. It is often studied phenomenologically with various QCD Monte- 
carlo models that have been tuned to fit the data 13] . Underlying event has been studied 
extensively at the Tevatron energy [T3 . Similar studies at the RHIC are very useful in 

understanding it's dependence on y/s, and can provide valuable constrains on the underlying 
event physics at the LHC. 

Figure lTTI shows the jet pair distribution in p + p collisions, reproduced from FigureOl The 
pedestal in the A(f> correlation, which represents the underlying event contribution, decreases 
quickly and becomes negligible at PT,assoc > 2 GeV/c. However, the level corresponding to 
minimum bias p + p events, denoted by the thick horizontal line, seems to decrease even faster. 
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Figure 10. Away side conditional yield as function of pr,trig for different ranges of xe in a) 
minimum-bias d + Au collisions and b) p + p collisions. 



Since minimum bias event has small hard-scattering contribution, the relative abundance of the 
pedestal in triggered events over the minimum bias events indicates that most of the underlying 
event comes from the initial or final state radiation of the hard-scattered partons. 
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Figure 11. Corrected condition yield in for p + p collisions (from FigureOJ). The thick 
solid line represents the average level for minimum bias events, i.e. it is equal to Yield pp /(27r). 



The underlying event at RHIC can be checked in QCD Monte-carlo models. We use the 
PYTHIA6.131, which seems to be able to reproduce the jet conditional yield as shown in 
Figure IT21 The roles of the initial/final state radiations are studied by switching them on and 



off in PYTHIA simuiation. Figure lT3k shows a typical ir^ — azimuthal correlation with (top 
histogram) and without (bottom histogram) radiation from the simulation. There is a significant 
enhancement in the pedestal level when radiations are enabled. We can perform a more 
quantitative comparison by plotting separately the jet contribution (double gauss component), 
the underlying event (the constant component) and the minimum bias event level as function 
of pt in Figure lT3b . The hierarchy of the three contributions can be clearly seen. For event 
tagged with a high jet, the spectra for both the jet fragmentation and the underlying events 
are much harder than that from typical minimum bias events. Current statistics from p+p does 
not allow a quantitative comparison with the models yet, a much larger dataset collected from 
recent RHIC p + p run in 2005 should help to address this question in the near future. 
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Figure 12. Conditional yield with trigger pion in 5 < prtrig < 6 GeV/c, compared with 
PYTHIA 6.131 for both same side (left panel) and away side (right panel). 
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Figure 13. PYTHIA MC simulation: a)A0 distribution for 3 < PT,assoc < 5 with (top 
histogram) and without (bottom histogram) radiation, b) the pt spectra for jet pairs and 
underlying event in p + p events with a > 5 GeV/c trigger, compared with minimum bias event 
yield. 



What about the underlying event in d + Au? FigureEl indicates that the underlying event 
levels are larger than those in p + p, although the properties of the jets are quite similar. Since 
the hard-scattering only happens in one of the nucleon-nucleon collision in d + Au, we can 



assume the ambient particle production mechanism is the same as in minimum bias nucleon- 
nucleon collision. In this case, the ambient particle production should simply scale as the nuclear 
modification factor, i?dAu measured in d + Au |ljj. Thus the underlying event yields in p + p 
and d + Au, UdAu arid U pp are connected to each other through the following simple relation, 

U dA u = U pp + R dAu (iV co u - 1) Yield pp (12) 

where Yield pp represents the hadron yield per event in minimum bias p + p collisions. Divide 
both side by Yield pp , we get, 

AdAu = A pp + -RdAu (-Wcoll - 1) (13) 

AdAu = t^dAu/Yieldpp, A pp = U pp /Yieldpp (14) 

note App denotes the ratio of underlying event yield to minimum bias event in p + p, which 
should be larger than 1 according to Figure lTTI and Figure IT^l 

4. Conclusion 

The di-jet decay kinematics are studied using tt^ — correlation in p+p and d + Au collisions. 
Measured jet widths, the calculated kx and distributions of p ou t are very similar between p + p 
and d + Au, which indicate no or small broadening in cold nuclear medium. Jet yield distribution 
in associated hadron px and xe are also similar between p + p and d + Au, consistent with no 
significant increase in jet multiplicity in d + Au relative to p + p. The dependence of the xe 
distribution on trigger px is weak in measured trigger px range. The underlying event yield in 
p + p is studied in PYTHIA Monte-carlo. Events containing a large px trigger appear to have 
a underlying event spectra much harder than the minimum bias hadron spectra. 
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